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The design of non-centrosymmetric structures usiaa 
copolymers has been a subject of great recent interest Jd'El 
Non-centrosymmetric bulk copolymer structures i-bave 
been analyzed using mean field continuous models .Eru In 
this communication we determine the conditions for ob- 
taining non-centrosymmetric structures, and how to tune 
the net macroscopic polarizibility solving a lattice model 
exactly. The model describes the self-organization, 
nanoaggregates of rodcoil molecules into polar films .B 
An important property of these nanoaggregates is their 
organization into macroscopic polar materials when cast 
from solution and annealed without the involvement of 
an external electric field. Transmission electron micro- 
graphs of film cross sections and x-ray scattering exper- 
iments reveal the formation ofJayered domains with a 
head-to-tail polar arrangement .H Interestingly, the mea- 
sured macroscopic polarization of these films is much less 
than one would-, expect for a monodomain of oriented 
nanoaggregates.B One possible explanation for the rel- 
atively small macroscopic polarity is cancellation among 
domains in the bulk of the film. In this communica- 
tion we explore the possibility of a non-randomly ori- 
ented microstructure by finding the ground state of a sim- 
ple model. The model is constructed to account for the 
competing interactions among nanoaggregates. The non- 
centrosymmetric structure of the nanoaggregates sug- 
gests they have a net dipole moment, and this leads to 
dipolar interactions among them. On the other hand, the 
enthalpic and entropic factors associated with contacts 
between coil and rod portions of neighboring aggregates 
suggests Ising-like nearest neighbor interactions. 

In our model the nanoaggregates are represented by 
dipoles of constant strength D on a cubic lattice of size 
a. While it is possible to consider very general models in 
which dipoles are free to orient in any direction, based 
on experiments the system we are interested in model- 
ing contains dipoles than can be aligned either parallel 
or anti-parallel to the z axis, perpendicular to the x — y 
substrate plane. Also, even though we included a penalty 
for reversing the orientation of a dipole with respect to 
its neighbors by means of an Ising coupling of strength 
J, there should be other steric forces associated with the 
shape and size of the nanoaggregates which are not con- 
sidered here explicitly. These forces could maintain the 
dipoles oriented along a given axis. The ground states of 
the dipole (J = 0) and of the Ising {D = 0) interactions 
along the z — x plane are shown in Figures l.a and l.b, 
respectively. The ground state when J = is antifer- 
romagnetic along x and y (columns along z), and when 

= is a homogeneous state with all dipoles pointing in 



the same direction (monodomain). With nonzero values 
of D and J, an intermediate stripe state of periodicity 
A is possible in which dipoles have the same orientation 
within equal size domains (A| = Aj^ = A/2) along the x 
or y axis and contiguous domains have anti-parallel ori- 
entations (see Figure Ic). This stripe structjire has been 
obtained pmiiiously in 2d magnetic systemaZI and in lipid 
monolayersBEl. We extended the model to finite thickness 
films .lij In infinite thickness films, we found a first order 
transition from a monodomain (A/a — oo) to anti-parallel 
domains with extremely small domain sizes (A/a < 4) as 
D/ J increases. However, in films of finite thickness, the 
domain periodicity A was found to decrease continuously 
as D / J increases. We add here a short range interaction 
between the dipoles and the surface and show that in this 
case the dipole up (A|) and dipole down (A|) domains no 
longer have equal widths as shown in Figure Ic, leading 
to net macroscopic polarization. We find that A and the 
macroscopic polarization can be modulated by varying 
the film thickness. 
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FIG. 1: Ground state configurations for a system of 
supramolecular aggregates with dipole-dipole interactions 
only a) , and with Ising interactions only b) . The stripe struc- 
ture under the surface influence c). 



Consider a 3d lattice, with dimensions La,, Ly and 
L^, composed of dipoles s^ with orientations up Sm — 
(0,0,1) or down s^ — (0,0,-1). The energy from the 
Ising interactions is given by. 



m ^="1X1 
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where < m, m' > are the nearest neighbors. The energy 



due to dipolar interactions is 
Ed 
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where f is the unit vector in the direction of m — m'. 

The monolayer (L^ = 1) case can be analyzed read- 
ily given that the second term in Ed vanishes. The 
existence of anti-parallel domains (stripes) in this 2d 
case has been justified by expanding the free energy 
in powers of the Fourier components (?!)(k) of a con- 
tinuous, field 4){v) proportional to the local polarizabil- 
ity asja F = J2k(t)i'^){Gi + G^)(f>{~k)+ quartic lo- 
cal terms. Here, the Ising contribution Gq is equal to 
Gq = 47J(— 1 -t- 2a?k'^) -\- ksT, where 7 is the number 
of nearest neighbors, T is temperature and /c^ is the 
Boltzman constant. The dipolar contribution is given 
by the 2d Fourier transform of the potential in eq 2, 
G^(fc) = (27ri?/a)(iF2(-l/2;l/2,l;-a2fcV4) - |fc|a), 
where if 2 is the generalized hypergeometric function. 
When D/J is small the — |fc| term in G^lka ^ 1) = 
(27rL)/a)(l - \k\a + k^a^ ji) added to the +k^ terms in 
Go gives a minimum F at a fc* mode that changes con- 
tinuously from fc* = (a monodomain) at DjJ — 
to fc* (a periodic structure) as socju-jas DjJ ^ 0. 
The lowest energy structure is a stripdld'Ej along x or y 
with periodicity A/a = 27r/fc*a = 47r + lQa/{D/J). This 
continuous analysis, however, is not accurate for neither 
D/J >nO-l nor low T.H The stripe structure melts via. 
defectsja with a local A given by the exact T = value,Ll 
A = Aexp(2aJ/D). The exponential decay of A/a from 
00 to 2 as D/J increases (followed by an antiferromag- 
netic state), agrees well with the exact numerical T = 
results given below. 

In order to analyze the stripe structure in films of ar- 
bitrary thickness we consider the system as a set oiy — z 
planes of one lattice thickness, stacked along x. Since 
within each plane all the dipoles have the same orien- 
tation, one can assign a new parameter to represent 
the configuration of all the dipoles within the plane. The 
Ising and dipolar energies per dipole between two planes 
for a repeat unit of width A is given by 
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respectively, where Vj^{\i' — i\) is the dipolar potential 
between planes i and i' . We compute V^(|«' — i\) by 
summing over all the individual dipole-dipole interac- 
tions that form the planes. 



(5) 



where mi and mj/ are the spins in planes i and i' respec- 
tively, has a strong dependence on L^, such that as 
increases Vq decreases at short distances and increases 
at long distances. Therefore the range of interaction is 
longer in thicker films. 

The interaction with the substrate can be easily in- 
cluded into this model. Let us assume the presence of 
a substrate in the x — y plane which favors the up con- 
figuration for the first layer of dipoles. The energy per 
dipole due to the interaction with this substrate can be 
written as. 
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This potential distorts the relative width of domains of 
dipoles up and down, such that A| /A 7^ 0.5. 
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FIG. 2: Results for A| vs. A for different surface interaction 
strengths for a film of 21 layers, S/J = (•), S/J=1(B) and 
S/J = 2(4). The inset shows A vs. D/J for the same S/J 
values. 

Figure 2 shows the effect of surface interaction on the 
ground state configuration, obtained by numerically min- 
imizing the total energy of the system. The calculations 
are carried out for a film of 21 layers of dipoles and infi- 
nite dimensions in the x and y directions. Surfaces with 
S/J = 0,1 and 2 give similar A variations as D/J in- 
creases. The domain size decreases continuously from 
infinite to finite values over a very narrow D/J range as 
shown in the inset in Figure 2. The effect of the surface 
interaction is revealed in the A| vs. A plot. The S/J = 
case yields A|/A ~ 0.5. As the surface interaction is 
turned on, the relative width of up and down domains 
changes such that A|/A > 0.5 (Figure Ic) and increase of 
the surface potential raises this difference further. 

The thickness of the film has a drastic effect on the 
ground state configuration, such that in the extreme limit 
of a bulk system the domain structure is destroyed. We 
analyzed three systems with film thickness, = 1 (a 



monolayer), = 21 and Lz = 201, all with S/J = 1. 
As shown in Figure 3, these three systems show a remark- 
ably different parameter dependence. For large values of 
D/J the monolayer has a periodicity of A/a = 2, and 
the domain size grows rapidly at about D/J^ 0.45, as 
shown in the inset in Figure 3. This sharp onset disap- 
pears for films with 21 and 201 layers. When we compare 
the ratio A|/A, it is clear that the monolayer acquires a 
monodomain configuration with full macroscopic polar- 
ization and the amount of polarization decreases as the 
film thickness increases. Since the model assumes that 
the surface interacts only with the first layer of dipoles, 
it is clear that the overall effect of the surface becomes 
less dominant as increases. 




FIG. 3: The effect of film thickness on the ratio A| /A for 
a monolayer (•), 21 layers (■) and 201 layers (♦). The inset 
shows A vs. D/J for the same film thicknesses. 



In our model we assume an Ising interaction which 
favors parallel configuration both along the z axis and 
within the x — y plane. Based on short range energetic 
considerations, one might conclude that the configura- 
tion along z would be anti-parallel, since this would en- 
able to some degree the phase separation of coil and rod 
blocks. However this simplification in our model alters 
the ground state of the system only when D / J is very 
small. That is, if we include an Ising interaction favor- 
ing an anti-parallel configuration along the z axis, the 
ground state is a monodomain of bilayers when D/ J <^ 1 . 
When D/J increases, however, the parallel configuration 
along z, favored by the dipolar interaction, dominates 
even when the Ising interaction favors an anti-parallel 
configuration along z. 

We conclude that the ground state of a system com- 
posed of dipolar supramolecular aggregates that interact 
with a surface is a periodic domain structure with a net 
macroscopic polarization. Furthermore, one can control 
the magnitude of this polarization through variations in 
the dimension of the system and the strength of the sur- 
face force constant. Monte Carlo simulations for some 
hexagonal lattices and smectic-A structures with direc- 
tional order also give stripe structures with similar A de- 
pendence on film thickness. 
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